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Abstract 



We consider two random matrix ensembles which are relevant for describing crit- 
ical spectral statistics in systems with multifractal eigenfunction statistics. One 
of them is the Gaussian non-invariant ensemble which eigenfunction statistics 
is multifractal, while the other is the invariant random matrix ensemble with a 
shallow, log-square confinement potential. We demonstrate a close correspon- 
dence between the spectral as well as eigenfuncton statistics of these random 
matrix ensembles and those of the random tight-binding Hamiltonian in the 
point of the Anderson localization transition in three dimensions. Finally we 
present a simple field theory in 1+1 dimensions which reproduces level statistics 
of both of these random matrix models and the classical Wigner-Dyson spec- 
tral statistics in the framework of the unified formalism of Luttinger liquid. We 
show that the (equal-time) density correlations in both random matrix models 
correspond to the finite-temperature density correlations of the Luttinger liquid. 
We also present a mechanism of the finite-temperature generation with break- 
ing the translational invariance in space by a metric with the event horizon, 
similar to the problem of Hawking radiation in the black holes. 
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1.1 Introduction 

It was known since the pioneer's work by F.Wegner [I] that the eigenfunctions 
ipi (r) of the random Schroedinger operator Htpi{r) = i^^j(r) at the mobility 
edge i?j = corresponding to the critical point of the Anderson localization 
transition, possess the property of multifractality. In particular, at E = E m 
the moments of the inverse participation ratio: 

Pn(E)=J2T / (M^)\ 2n m-E))ocL- d ^- 1 \ (1.1.1) 

r i 

scale as a certain power-law with the total size L of the system. This power is 
a critical exponent which depends only on the basic symmetry (see Ref.[2] and 
Chapter 3 of this book) of the Hamiltonian H and on the dimensionality of space 
d. The true extended states of the Schroedinger operator are characterized by 
all d n = d. This allows to interpret d n < d as a certain fractal dimension which 
depends on the order n of the moment. As a matter of fact the statistics of 
critical eigenfunctions are described by a set of fractal dimensions d n which 
justifies the notion of mufti- fractality. 

Another aspect of criticality is the scaling with the energy difference lo = 
\Ei — Ej\ between two eigenvalues. It is similar to the dynamical scaling and is 
relevant for the correlation functions of different eigenfunctions ipi(r) and ijjj( r )- 
The most important of them is the local density of states correlation function: 

cv, r)=J2 Ed^wi 2 + R )i 2 S ( E - E ^ E + w - E i))- ( LL2 ) 

r ¥i 

This correlation function is relevant for the matrix elements of the two-body in- 
teraction. The dynamical scaling connects the power law behavior of C(u>, 0) oc 
lo^^ with that of C(0, R) oc R~( d ~ d2 ) by a conjecture [3] on the dynamical 
exponent: 

R d ^u;, n = l-!j. (1.1.3) 

Although there is an extensive numerical evidence in favor of conjecture Eq. (|1.1.3j) 
its rigorous proof has been lacking so far. 

Last but not least, there is a growing interest to the spectral (level) statistics 
in quantum systems whose classical counterparts are in between of chaos and 
integrability [3], the simplest of them being the two-level spectral correlation 
function (TLSCF) R{uj) oc ^] R C(w,fl). Some of them apparently share the 
characteristic features of spectral statistics at the Anderson localization tran- 
sition, e.g. the finite level compressibility < x < 1 [SHU], and the Poisson 
tail of the level spacing distribution function P(s) oc e~ s / 2x combined with 
the Wigner-Dyson level repulsion P{s) oc s^ at small level separation s <C 1 
[7]. There is a conjecture that these features are also related with the multi- 
fractality of critical eigenfunctions, however its exact formulation has not been 
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developed beyond the limit of weak multifractality x 1, where one can show 
that x = \{l-d 2 /d) [31]. 

The main reason for scarcity of rigorous knowledge about multifractality of 
critical eigenfunctions is a lack of exactly solvable models of sufficient generality. 
The most popular three-dimensional (3D) Anderson model of localization [8] is 
quite efficient for numerical simulations but so far evaded any rigorous analytical 
treatment in the critical region. More perspective seemed to be the Chalker's 
network model for the Quantum Hall transition [ L 9j and its generalizations. 
However, numerous proposals for the critical field theory were not successful 
so far [10]. In this situation the random matrix theory may prove to be the 
simplest, universal and representative tool to obtain a rigorous knowledge about 
the critical eigenfunction and spectral statistics 

1.2 Non-invariant Gaussian random matrix theory 
with multifractal eigenvectors. 

Guided by the idea that multifractality of eigenstates is the hallmark of critical- 
ity, we introduce the Gaussian random matrix ensemble [12} [TT] which eigen- 
vectors obey Eq. tjl.l.ip with L being replaced by the matrix size N. This 
random matrix theory and its modifications describes very well not only the 
critical eigenfunction statistics at the Anderson localization transition in three- 
dimensional (3D) Anderson model but also the off-critical states close to the 
transition [13]. The critical random matrix ensemble (CRMT) suggested in 
|12[ [TT] is manifest non-invariant, and is defined as follows: 



where H nm is the Hermitean N x N random matrix with entries H nm (n > 
m) being independent Gaussian random variables; (3 = 1,2,4 for the Dyson 
orthogonal, unitary and symplectic symmetry classes, and b > is the control 
parameter. This CRMT can be considered as a particular deformation of the 
Wigner-Dyson RMT which corresponds to b = oo. 

As is clear from the definition Eq. (jl.2.1j) the variance (|^nm| ) is non- 
invariant under unitary transformation H — ► U HU^ . The existence of the 
preferential basis is natural, as this RMT mimics the properties of the An- 
derson model of localization which happens in the co-ordinate space, and not 
necessarily e.g. in the momentum space. 

The critical nature of the CRMT is encoded in the power-law decay of 
the variance matrix, the typical off-diagonal entry being proportional to \n — 
m\~ l in the absolute value. This is exactly the decay with the power equal 
to the dimensionality of space (d = 1 in the case of matrices). In contrast to 




(1.2.1) 
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the Wigner-Dyson RMT which probability distribution is parameter-free, the 
CRMT is a one-parameter family. The parameter b controls the spectrum of 
fractal dimensions d n = d n (b). One can show [12\ [T4] that both at b 3> 1 and 
b -C 1 the scaling relation Eq. (ll.l.i~j) holds true, and the basic fractal dimension 
is equal to: 

\ cpB, S<1 17 

4-1 1 /3 1 

where cp = - 2 4 , and B = b-K^ 2i. 

Eq. (jl.2.2p can be cast in a form of the duality relationship: 

d 2 {B) + d 2 {B- 1 ) = 1. (1.2.3) 

The duality relation has been checked numerically [15] for the unitary CRMT 



1.10 




B 



Figure 1.1: Numerical verification of the duality relation Eq. (jl.2.3p . done by 
the box counting method [15] for the unitary CRMT ensemble Eq. (ll.2.ip . The 
deviations from 1 in Eq. (jl.2.3p do not exceed 1% in the entire region of B G 
[0,1]. 

Eq. (|1.2.ip . The results are presented in Fig. 1.1. In particular, it was found 
that 2d 2 (B = 1) = 1.003 ± 0.004. The fact that the deviation of the sum 
d 2 (B) + d 2 {l/B) from 1 does not exceed 1% in the entire region of B G [0, 1] 
looks extremely interesting. However, it is not known yet whether an exact 
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function d,2(B) obeys this remarkable duality relationship which suggests that 
d,2 = h at b = —r^ — . 

There was vast numerical evidence that the CRMT Eq, (|1.2.ip reproduces 
the main qualitative features of the critical eigenfunction and spectral statistics, 
including the power-law behavior of the DoS correlation function C(u, 0) [see 
Ref.|13] and references therein], distributions of moments of inverse participa- 
tion ratio P q [14], role of rare realizations and multifractality spectrum close to 
termination point [16], and the hybrid Wigner-Dyson & Poisson level spacing 
distribution [7]. 

The very fact that by a choice of control parameter B one can fit quite 
accurately the critical statistics of both eigenvalues and eigenfunctions of the 
Anderson localization model, is extremely encouraging. 



1.3 Invariant RMT with log-square confinement 

Quite remarkably, there is an invariant (but non-Gaussian) RMT whose two- 
level spectral correlation function Rn(s,s') is closely related with that of the 
CRMT discussed in the previous section. Namely, in the unfolded energy vari- 
ables s in which the mean spectral density ("global Density of States") is equal 
to unity, and in the large B limit one finds: 

Boo{s S ) |inv. — Roo(s S ) |non— inv. j (1.3.1) 

where i?oo(s - s') = lim s ^ 00 (limAr^ 00 Rn(s, s')). 

For the unitary symmetry class one can show [11] that: 

Tj f l\\ 1 2 2 sin 2 (7T(g-g / )) (3 

Roo{S ~ S ) |non-inv = 1 ~ IT K ■ — , K =- X (b), (1.3.2) 

sinn [7t z k(s — s'JJ ^ 

where k is related with the level compressibility x(b) (see Eq. (jl.8.3[) . (jl.8.6|) 
below). This suggests [2] the form of the kernel K(s — s') of the invariant 
RMT which is the only input one needs to compute all many-point spectral 
correlation functions: 

K ( . .y^'V f 1 - 3 - 3 ) 

smn[7r z K(s — s )\ 

Thus a remarkable correspondence between an invariant and a non-invariant 
RMT can be conjectured [TT] which allows to use the full power of the unitary 
invariance for calculation of spectral statistics. 

Now it is time to specify the invariant RMT whose counterpart in Eq. (jl.3.ip 
is the critical RMT with multifractal eigenstates. The probability distribution 
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for the random matrix Hamiltonian H is |17l [TT] : 

oo 

P(i?) ocexp -fitrV(H) , V{x) = ^ In [l + 2q n+1 (1 + 2x 2 ) + g 2(n+1) 

n=0 

(1.3.4) 

where 1 < g < is a control parameter. 

It is extremely important that the "confinement potential" V(x) at large 
|x| behaves like: 

V(x) -> A In 2 ^ = — — — . (1.3.5) 

ln(g x ) 

thus being an example of a "shallow" confinement. The correspondence Eq. (jl.3.1j! 
holds for lng -1 << 1, while for lng -1 > tt 2 an interesting phenomenon of "crys- 
tallization" of eigenvalues happens with the TLCF taking a triangular form 

The form of the confinement potential in Eq. (|1.3.4p is quite specific, even 
within the class of shallow potentials with a double-logarithmic asymptotic 
behavior Eq. (jl.3.5p . It has been chosen in Ref.|17j in order to enable an exact 
solution in terms of the g-deformed Hermite polynomials. However, there exists 
another, more simple argument why this particular form leads to an exact 
solution |19j . The point is that the measure Eq. (jl.3.4p can be considered as a 
generalized Cauchy distribution: 

oo N n+1 

v(Ei) = TT —-- "TT — — r n = — \. . (1.3.6) 



1.4 Self- unfolding and not self- unfolding invariant 
RMT. 

As we will see below, the cause of all the peculiarities of the RMT with double- 
logarithmic confinement is the fact that the mean density of states p oc (£ ; ) = 
lmi/v^oo ^2 i= i(S(E — Ei)) approaches some stable non-trivial form in the bulk of 
the spectrum as the size of matrix N tends to infinity. This is in contrast to the 

Wi gncr-Dyson classics! RMT where Hni^v >oc Pn(E)/V2N = 1 is independent 

of E. This is the reason why the Wigner-Dyson RMT can be referred to as 
self-unfolding, while the RMT with the double-logarithmic confinement is not 
self-unfolding. 

In order to find a criterion for an invariant RMT to be self-unfolding, let us 
apply the Wigner-Dyson plasma analogy [2]. According to this approximation 
the mean density of states pn(E) obeys the integral equation of the equilib- 
rium classical plasma with logarithmic interaction subject to the confining force 
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-dV/dE: 

dE' dV 



v.p.j PN (E') ¥ - w = - = f(E), (1.4.1) 



oc 



where v.p. denotes the principle value of the integral. The solution to this 
equation in the case of symmetric confining potential is: 

Pn(E) = -i y/D% - E* v.p. f{E ' ] (1.4.2) 

where the band-edge Dn should be determined from the normalization condi- 
tionf°g N p N (E')dE' = N. 

One can see that if \f(E')\ increases slower than \E'\ as \E'\ — * oo, the 
integral converges to a non-trivial function as N and Dn — > oo: 

Poo(E) = \ v.p. [ + °° dE' Jp*L. (1.4.3) 
" J-oo E' - E 

Otherwise, the large values of \E'\ ~ Djy dominate the integral in Eq. (jl.4.2p . 
so that at a fixed E and N — > oo (when Djy ^> \E\) the dependence of Pn(E) 
on E disappears. We conclude therefore [20^ [2T] that the criterion of a non- 
self-unfolding RMT is: 

lim ^M = 0. (1.4.4) 

\x\— »oo \x\ 

Let us consider the shallow confining potential with the power-law large- x 
asymptotic behavior 

V(x) = A\x\ a , (a<l). (1.4.5) 
Then at large \E\ > 1 we have [20||2"T]: 

Act /ira\ 1 . , , 

Poo{E ) = — t & n [—) (1.4.6) 

In particular for the log-square confining potential 

V{x)= A\n 2 \x\= AXvn ^ =^lim^ ' 1 — . (1.4.7) 



one finds using the linear dependence of Eq. (jl.4.3p p^ (E) on f(E) and Eq. (jl.4.5p : 

Poo(E) = A. (1.4.8) 

There is a qualitative and far-reaching difference between the shallow power- 
law confinement potential with < a < 1, and the log-square confinement. 



""^For V(x) oc ln d |k| with d > 2 the leading term at large _B is pao(E) oc ln E 



In""" 8 \E\ 
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Although both lead to a non-self-unfolding RMT, the case of finite a can be 
called weakly non-self unfolding, because for large enough E the variation of 
the mean density at a scale of the mean level spacing A = p^ 1 is much smaller 
than the density itself: 

Poo (E + A) - Poo (E) _Pk (xlEr ^ (1A9) 



Poo{E) p, 



OC 



In contrast to that, the log-square confinement is strongly non-self unfolding, 
since the ratio in Eq. (|1.4.9j) is always finite H. 

We conclude that depending on the steepness of the confinement potential 
at large E there are three qualitatively different cases [201 EI] : 

• self- unfolding RMT for lim^^oo V(x)/\x\ > 

• weakly non- self- unfolding RMTioi lim^^oo V(a;)/|s| = but 3a > such 
that lim^ooV^xVlxl" > 

• strongly non-self unfolding RMT if Va > holds lirn r _+ 00 l/(x)/|x| a = 

The first case is characterized by the Wigner-Dyson universality of the spectral 
correlation functions. In the second case this universality holds only approxi- 
mately for sufficiently large distance from the origin, while if one or two energies 
are close to the origin, the Wigner-Dyson universality is no longer valid, even 
after unfolding. In the third case, the Wigner-Dyson universality is not valid 
also in the bulk of the spectrum. 



1.5 Unfolding and the spectral correlations 

Let us consider the power-law confinement Eq. (|1.4.5p and find the correspond- 
ing unfolding co-ordinates s(E) in which the spectral density is 1: 

s(E) = sgn(E) Jj Poo(E') dE' = - tan (^) sgn(£) \E\ a . (1.5.1) 

Note that for large enough E the unfolding co-ordinates are given by Eq. (11.5.I]) 
even if V(E) is not a pure power-law at small E but is rather deformed to have 
a regular behavior at the origin. 

The corresponding unfolding co-ordinates for the log-square potential Eq. (|1.4.7|) 

are: 

s(E) = A sgn(£) ln(c\E\), cE(s) = sgn(s) e^. (1.5.2) 

where c is a constant which depends on the regularization of the log-square 
potential close to the origin. 

2 We put the confinement potentials ln d \E\ with d > 2 to the class of strongly non-self 
unfolding potentials even though the ratio Eq. (|1.4.9|l logarithmically vanishes at large E. 
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As we will see, the exponential change of co-ordinates Eq, (jl.5.2p leads to 
dramatic consequences for two-level spectral correlations, and has even a far- 
reaching analogy in the physics of black holes [2Q|, [22] . 

In order to show how a non-trivial unfolding changes the form of the spec- 
tral correlations consider the spectral kernel Kn(E, E') given in terms of the 
orthogonal polynomials by the Christoffel-Darboux formula [2]: 



K N (E,E') 



(p N - X {E)(p N {E') - y N - X {E')ip N (E) 



E-E' 



;i.5.3) 



where ipN(E) is the (properly normalized) " wave- function" which is related 
with the orthogonal polynomials p n : 

r+oo 

Vn{E) = Vn {E) e~^l\ \ p n {E')p m (E') e^ v ^ dE' = S nm . 



(1.5.4) 

Any spectral correlation function can be expressed [2] in terms of the kernel 
Eq. (jl.5.3p . In particular, the n-point Density of States correlation function at 
= 2 takes the form [2]: 

n 

R 00 (E 1 , ...E n ) = det {K^E^Ej)} JJtf- 1 (£;,£;). (1.5.5) 

8=1 

Below we derive a semi- classical form of the kernel which is valid in the N — > oo 
limit provided that the coefficient A in Eqs. (jl.4.5p . (jl.4.7p is large. In particular 
for the RMT with log-square confinement given by Eq. (jl.3.4p the condition of 
applicability of the analysis done below is lng -1 < 2ir |17j . 

In this semiclassical limit the "wave functions" take the form (we assume 
that the confinement potential is an even function of E) : 

ipN-i(E) = sm(7Ts(E)), <fN = cos(tts(E)), (1.5.6) 

if N is even and cos — * sin if N is odd. Then one finds in the unfolding 
co-ordinates: 

Equation Eq. ()1.5.7p demonstrates that as long as the semiclassical approach 
applies, the non-trivial unfolding E = E(s) is the only source of deformation 
of the kernel and its deviation from the universal Wigner-Dyson form. 

The main conclusion one can draw from Eq. (|1.5.7p is that the translational 
invariance is lost in the N — > oo limit for all non-self- unfolding invariant RMT. 
It is not sufficient to have the mean density poo( s ) equal to unity for all values 
of s in order to have all the correlation functions of the universal Wigner-Dyson 
form. 
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1.6 Ghost correlation dip in RMT and Hawking ra- 
diation. 

The break-down of translational invariance takes especially dramatic form in 
case of the log-square confinement. We will show below that in this case a ghost 
correlation dip appears in the two-level spectral correlation function which po- 
sition at s ~ — s' is mirror reflected relative to the position of the translational- 
invariant correlation dip at s ~ s'. 

Indeed let us consider the semiclassical kernel Eq. (jl.5,7p with the expo- 
nential unfolding Eq. (jl.5,2p for s>l and \s — s'\ « \s\. Then after a simple 
algebra we obtain for s s' > the Two-level Spectral Correlation Function given 
by Eq. (fTX2l) ((3 = 2): 

„ . , . 99 sin 2 (7r(s — s')) 1 lno -1 

sinrr [7t 2 k(s - s')] 27T A 4ir z 

(1.6.1) 

This correlation functions exhibits a dip of anti-correlation (level repulsion) 
at small s — s' and approaches exponentially fast the uncorrelated asymptotic 
Roo(s, s') = 1 for \s — s'\ > 1/tt 2 k. An amazing fact is that the anti-correlation 
revives again when s + s' becomes small. Considering s s' < and plugging the 
exponential unfolding Eq. (jl.5.2p into Eq. (jl.5.7p we obtain for the unitary case 
(3 = 2: 

„ , , k 99 sin 2 (7r(s — s')) 

Roc{ss' <0) = 1-7T 2 K 2 i2 ; 9 V ^— . (1.6.2) 



This is the ghost anti- correlation peak discovered in Ref.[20j but also present in 
the exact solution [T7] . 

The existence of such a ghost correlation dip is requested by the normal- 
ization sum rule [20] which for the invariant RMT survives taking the limit 
N -> oo: 

/+oo 
(1 - Roo^s, s'))ds' = 1. (1.6.3) 
-oo 

Substituting in Eq. (jl.6.3p the sum of the translational invariant and the ghost 
peak found from Eqs. (|1.6.1jl . (|1.6.2j> and doing the integral one obtains: 

coth(/t _1 ) — — — cos(4vr s) w 1. (1.6.4) 

smh(K 1 ) 

This expression is equal to 1 with the exponential accuracy in e" <C 1, which 
is exactly the accuracy of the semiclassical approximation Eqs. (|1.6.ip . (|1.6.2p . 

It appears that the situation with the ghost correlation dip has an important 
physical realization [22]. Let us consider the sonic analogue of the black hole 
[23\ I24j described in Fig. 1.3. The stream of cold atoms moves with the constant 
velocity vq, while the interaction between atoms is tuned so that the sound 



1 . 6. GHOST CORRELATION DIP IN RMT AND HAWKING RADIATION. 1 1 




Figure 1.2: The irreducible part Roo(s,s') - 1 of the TLSCF obtained [21] 
by the classical Monte-Carlo simulations at a temperature T = i on the one- 
dimensional plasma with logarithmic interaction in the log-squared confinement 
potential Eq. (|1.4.7|) with A = 0.5. Two dips of the anti-correlations are well 
seen with the smaller ghost dip shown in detail in the insert. The edges of the 
correlation function correspond to the finite number of particles (energy levels) 
N = 101. 



velocity is larger than the stream velocity for x < —5 and smaller than the 
stream velocity for x > 5. The latter region where no phonon may escape from, 
is analogous to the interior of a black hole. The point in space where c(x) = vq 
is analogous to the horizon in the black hole physics and the phonon radiation 
emerging when the horizon is formed, is analogous to the Hawking radiation 
[25]. 

As it was shown in Ref.[24j, the Hawking radiation possesses a peculiar 
correlation of photons (phonons): despite being bosons, they statistically repel 
each other forming a dip in the density correlation function not only at x ~ x' 
but also near the mirror point x ~ —x', with the envelopes of the normal and 
ghost dips proportional to sinh~ 2 and cosh -2 , respectively. As in the case of 
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horizon 
region 



Super-sonic flow V >c 



Exterior of "black hole 



- o 



Interior of "black hole" 



Figure 1.3: The sonic analogue of a black hole (courtesy of I.Carusotto): the 
stream of cold atoms with the the velocity vq and the phonon velocity c(x) 
changing near the origin from the value c\ > vq to the value C2 < Vq. This 
region of the super-sonic flow is analogous to the interior of the black hole 
where light cannot escape from. 

RMT with log-square confinement, the ghost correlation dip appears in the 
absence of any symmetry that would require the eigenvalues E{ or the Hawking 
bosons to appear in pairs simultaneously at points rbJSj (±£j). More close 
inspection shows that the mechanism of its appearance in quantum gravity 
and in the random matrix theory is very similar: it is the exponential change 
of variables similar to Eq. (|1.5.2|) which is aimed to make flat the mean DoS 
(unfolding) in RMT or to make flat the metric of space-time in the quantum 
gravity (exponential red-shift) [22]. Thus formation of the strong non-self- 
unfolding RMT with the log-square confinement appears to be analogous to 
formation of a horizon in the general relativity. 

1.7 Invariant-noninvariant correspondence. 

Note that the translational-invariant part of spectral correlations described by 
Eq. (jl.6.ip is well separated from the translational non-invariant ghost dip if 
|s| S> |s — s'\. As correlations exponentially decrease at |s — s'\ > (7t 2 k) _1 , 
the scale separation takes place when \s\ (7t 2 k) _1 . Essentially what happens 
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because of the scale separation is that (like in a black hole) the world is divided 
in two parts, and in each of them local spectral correlations (which are approx- 
imately translational- invariant ) are not affected by the presence of a "parallel 
world". Moreover, the local correlations in the invariant RMT with log-square 
confinement and k <C 1 were conjectured (for the two-level correlations this 
conjecture has been proven in [11]) to be the same as in the non-invariant crit- 
ical RMT |11^ 126]. This conjecture is related with the idea of the spontaneous 
breaking of unitary invariance in RMT with shallow confinement [20] . Although 
this idea has never been proven or even convincingly demonstrated numerically, 
it seems to be the only physically reasonable cause of the invariant-noninvariant 
correspondence like the one given by Eq. (|1.3.ip , 

The invariant-noninvariant correspondence allows to use the invariant RMT 
with log-square confinement for computing spectral correlations of non-invariant 
RMT considered in sec. 1.2. The idea of such calculations is to use the kernel 
in the form Eq. (jl.3.3p to be plugged into the conventional machinery of invari- 
ant RMT [2j. Here we present some results of such calculations taken from 
Rcf. [27]. One can see from Fig. 1.4 that the choice of one single parameter A in 




Figure 1.4: Level spacing distribution function (after Nishigaki |27| ) for the 
orthogonal and the unitary ensembles with log-square confinement (black solid 
lines) and the corresponding distributions for the 3D Anderson localization 
model with critical disorder (data points). Grey solid lines are the Wigner- 
Dyson and the Poisson distributions. The same values of A = 0.34 (orthogonal 
ensemble) and A = 0.14 (unitary ensemble) allow to fit well both the body and 
the far tail of the distribution shown in inserts. 

the invariant RMT with log-square confinement allows to fit well both the body 
and the exponential tail of the level spacing distribution function obtained by 
numerical diagonalization of the 3D Anderson localization model at criticality. 
The corresponding parameters b of the non-invariant CRMT Eq. (|1.2.1j) should 
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be chosen so that the parameter k in the kernel Eq. (jl.3.3p expressed through 
A and b is the same. Comparing expressions for k in Eqs. (|1.3.2|) . (jl.6.ip we 
conclude that b is a solution to the equation: 

' Px(b). (1.7.1) 



7T 2 A 



For the unitary ensemble the function x(b) in the range x = (2vr 2 yl) _1 ~ 0.36 is 
well described [M] by the large-6 asymptotic formula [14J x(6) = (47r6) _1 . This 
gives an estimate b = 0.22. Numerical diagonalization of the CRMT shows that 
b = 0.22 corresponds to d<i ~ 0.46 ±0.01. To establish the correspondence with 
the Anderson localization model of unitary symmetry at critical disorder we 
note that the fractal dimension in the Id CRMT corresponds to the reduced 
dimension dij'S of the three-dimensional Anderson model. Thus the fractal 
dimension of the Anderson model should be compared with 0.46 x 3 = 1.38. 
It appears to be in an excellent agreement with the direct diagonalization of 
the Anderson model of unitary symmetry which gives [35] c?2 = 1-37. This 
example demonstrates that the values of parameters A and b found from fitting 
the spectral statistics of RMT to that of the 3D Anderson localization model, 
automatically give an excellent fitting of the eigenfunction statistics. 

1.8 Normalization anomaly, Luttinger liquid anal- 
ogy and the Hawking temperature. 

As has been already mentioned, the local spectral correlations with \s — s'\ <C \s\ 
are well described by the translational- invariant kernel Eq. (|1.3.3p . An im- 
portant difference between this kernel and that of the conventional Wigner- 
Dyson theory is that it contains a second energy scale k _1 in addition to the 
mean level spacing A = 1. Moreover, the way this scale appears through s — 
s' — > sinh(7r 2 K (s — s')) suggests an analogy with the system of one-dimensional 
fermions at a finite temperature. Indeed, the density correlations of non-interacting 
one-dimensional fermions (whose ground state correlations are equivalent to the 
Wigner-Dyson spectral statistics at (3 = 2) is described by Eq. (|l,3.2p . where 

T = 7TK. (1.8.1) 

Given the kernel Eq, (|1.3.3p it is not difficult to obtain the level density corre- 
lation functions for f3 = 1,4 using the standard formulae of the Wigner-Dyson 
RMT [2j. It appears that both in the region of \s > 1 and in the region of 
\s — s'\ -C k^ 1 they coincide with the generalization of Eq. (jl.3.2p to /3 = 1, 4 ob- 
tained [26j directly from the non-invariant CRMT Eq. (jl.2.ip . For k <C 1 these 
regions overlap on a parametrical large interval, and one can again demonstrate 
the invariant-noninvariant correspondence. 
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It is important to note that breaking the unitary invariance explicitly by a 
finite b in Eq, (|1.2.ip leads to a normalization anomaly. Namely, the normaliza- 
tion sum rule [20] which exactly holds for the finite size of matrix N: 

r+oo 

2 (l-R N (s' -s))ds' = 1, (1.8.2) 

Js>0 

gets violated if the limit N — > oo is taken prior to doing the integral. For 
instance at (3 = 2 we have with exponential accuracy e~ x l K <C 1: 



V 



poo 

2 (1 - Roods' - a)) ds' = cothfc -1 ) - k « 1 - k. (1.8.3) 

Js>0 



>s>0 

For (3 = 1 the corresponding expression is: 

r) = 2 cotfa -1 ) - tanh(l/2«) - 2k « 1 - 2k. (1.8.4) 
This leads to the finite /ewe/ compressibility [51] : 

X = -tt ((n - n) 2 ) = 1 - t/, iV»n=(n)>l, (1.8.5) 
an 

and the exponential (instead of the Gaussian in the Wigner-Dyson RMT) tail 
of the level spacing distribution function: 

lnP(*)«-^-. (1.8.6) 

Both properties Eq. (|1.8.5p and Eq. ([1.8.6|) are the signatures of criticality [7]. 
Due to invariant-noninvariant correspondence the normalization anomaly Eq. (|1.8.3p 
holds also for the invariant RMT with log-square confinement thus again rasing 
a question on the spontaneous breaking of unitary invariance. 

In the absence of the formal proof of this conjecture, we present here a sim- 
ple theory which may unify both of the ensembles. The idea of such a theory 
stems from the fact [28J that the classical Wigner-Dyson RMT [2] is equiv- 
alent to the ground state of the one-dimensional Calogero-Sutherland model 
of fermions with inverse square interaction of the strength 5 ( £ — 1 ) in 



2 

a harmonic confinement potential. The large-scale properties of such a model 
are described by the Luttinger liquid phenomenology |3L)j . Its simplest finite- 
temperature formulation [32j is in terms of the fvee-bosonic field in a 1 + 1 
space-time <E»(s,r) = \ [$r(s,t) + $>l(s,t)] with the quadratic action: 

o f l/T f +<x> 

5 r[*] = r^/ dr ds[(d s <S>) 2 + (d T <S>) 2 ], $(s,t) = $(s,t + 1/T). 

(1.8.7) 
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The density correlation functions are given by the functional averages of the 
density operator: 

1 

p(s,r) = p + -d s $(s,T) + A 1 cos[2tts + 2^(s,t)]+A 2 cos[4vr s + 4$(s, r)] + 

(1.8.8) 

where A^ are structural constants which are determined by details of interaction 
at small distances and take some fixed values for the Calogero-Sutherland model 
corresponding to the symmetry class (3. Using Eq. (ll,8.8j) one may express the 
density correlation functions in terms of the fundamental Green's function: 

($( S , r) $( S ', t)) s ~ \(*{8, r) *( s , r)) s ~ r)) *{st , t)) s = ^ G(s, s'), 

(1.8.9) 

where (■■■)s is the functional average with the action S[$>]. For example, the 
two-level spectral correlation functions are equal to: 

Rvois, s') = 1 + \ d s d s , G(s, s') + ^—s- cos(2vr s) e 8 ^' 1 G(s ' s,) . (1.8.10) 

71 P /3(2tt 2 )^ 

For T = the Green's function is G(s — s') = — ln(s — s') 2 at \s — s'\ 1 
and one reproduces the asymptotic form of the Wigner-Dyson correlations. 

In order to reproduce by the same token the corresponding correlation func- 
tions [26] for the non-invariant CRMT Eq. (jl.2.ip one merely substitutes in 
Eq. (ll.8.7p the finite T given by Rq. (|1.8.ip and replaces the zero-temperature 
Green's function by the finite-temperature one: 

G(s - s') = — In ( — — — V (1.8.11) 

v 1 2vr V sinh ( 7rT l' s - s \)J 

Thus we conclude that the deformation of the Wigner-Dyson RMT given 
by Eq. (jl.2.ip with large but finite b, retains the analogy with the Calogero- 
Sutherland model. However, instead of the ground state, the non-invariant 
CRMT with k 1 is equivalent to the Calogero-Sutherland model at a small 
but finite temperature Eq. (|1.8.ip . 

It is remarkable that there exists a deformation of the free-boson functional 
Eq. (ll.8.7p capable of reproducing the two-level correlation function for the in- 
variant RMT with log-square confinement, including the ghost correlation dip. 
All one has to do for that is to replace the action Eq. (jl.8.7p defined on a cylinder 
by that defined on a curved space-time with a horizon: 

/■oo /*+oo P 

\ dr ds[(d s <S>) 2 + (d T $) 2 ] -» / #Zy/=rtg)grdp*d v $, (1.8.12) 

JO Joo J 

where g = det g^ u with g^ u being the metric, i.e. ds 2 = g^ v dt^ d£ u with 
x 1 = x and x° = t = —it. The main requirement for the metric is that 
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the transformation of co-ordinates to the frame (x, t) where the metric is flat 
("Minkovski space" with ds 2 = dx 2 — dt 2 ), is exponential for large enough x 
with a factorized x and t dependence, e.g. 

x = <p(x) cosh(t/^4), i = tp(x) smh(t/A), (1.8.13) 
ip(x) ~ sgn(x) at \x\ > A = 2vr T. (1.8.14) 

In order to represent the invariant RMT with an even confinement potential, 
the function (p(x) should be odd. This last requirement plus the continuity of 
the function at x = necessarily implies goo(x = 0) = 0, that is x = is the 
horizon. One example of such a metric is: 

x = A sinh(x/A) cosh(i/A), t = A swh(x/A) siah(t/A) (1.8.15) 
ds 2 = cosh 2 (x/,4) dx 2 - sinh 2 (x/,4) dt 2 (1.8.16) 

= | tanh(s/^)| {d x §) 2 + | coth(x/A)| (d r $) 2 . (1.8.17) 

It maps the strip x>0, Q < t < 2ir A and the strip x<0, < t < 2ir A 
separately onto the entire plane (x,f). 

In order to compute the Green's function on the curved space-time corre- 
sponding to Eq. (jl.8.13p we use the well known formula |33j : 

G(z, A = ~^ In \z(z) - z(z')\ - i- ln[\d z z(z) d z ,z(z')\], (1.8.18) 

where z = x ± if , z = x ± it and d z = d x =F id T with ± = sgn(x). Then from 
Eq. (jl.8.13p for x, x' sufficiently far from the origin one easily obtains: 



G(x,x') = -—In 



Mx)-y(xQ) s 
4\( P (x)(p(x')\ 



^J_(ln[2Asmh[{x-x')/2A]], x x' > . . 
2tt \ln[2^cosh[(x + x')/2^]] , xx' <0 ^ j 

The origin of sinh and cosh in Eq. (|1.8.19p is very much the same as that in 
Eqs. (jl.6.ip . (|1.6.2p . The exponential transformation of coordinates Eq. (jl.8.13|) 
plays the same role as the exponential unfolding Eq. (jl.5.2p . Finally substituting 
Eq. (jl.8.19p into the expression for the two- level correlation function Eq. (jl.8.10l) 
one reproduces Eqs. ()1.6.ip . (ll.6.2p . 

Note that the new scale A = (2ttT)~ 1 sets the temperature scale T given 
by Eq. (ll.8.ip which has a meaning of the Hawking temperature in the black 
hole analogy. While the finite temperature arises because of the periodicity of 
transformations Eq. (jl.8.15p in the imaginary time r = it, this compactification 
is different from the standard one: each of the semi-strips x > and x < 
are mapped on the entire plane (x, f) independently of each other. This is 
essentially the effect of the horizon. 
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1.9 Conclusion 

I would like to close this Chapter by some concluding remarks. First of all, 
it is by no means a comprehensive review but rather some introduction to the 
subject written by a physicist motivated by the physics applications and not 
by the formal rigor. My goal was to show that the subject is rich and poorly 
explored and that the efforts are likely to be rewarded by non-trivial discoveries. 
Let me formulate in the end the (highly subjective) list of open problems as I 
see them. 

• further study of the non-invariant Gaussian critical RMT 

Some progress has been made in this direction by development of the 
regular expansion in b <C 1, the so called virial expansion method [36] 137} 
[38] , which extends ideas of Refs.[39|[T4]. This approach made it possible to 
compute analytically the level compressibility up to b 2 terms and find an 
extremely good and simple approximation |38} [13] to the LDoS correlation 
function Eq. (|1.1.2j) . These works are the basis for a perturbative proof of 
the dynamical scaling Eq. (|1.1.3j) [5] and the duality relation Eq. (|1.2.3p . 

• N on- perturbative solution to the non-invariant critical RMT 

Some nice relationships such as the duality relation Eq. (|1.2.3p . raise a 
question about a possibility of exact solution to the CRMT. In our opinion 
this possibility does exist. 

• Invariant-noninvariant correspondence 

This is a very interesting issue with lots of applications in case the origin 
of this correspondence is understood. It is also important to invest some 
efforts to study the issue of the correspondence between the spectral and 
eigenvector statistics, in particular the possible spontaneous breaking of 
unitary invariance and emergence of a preferential basis [20 1. 

• Level statistics in weakly non-self-unfolding RMT 

This is a broad class of invariant RMT where the Wigner-Dyson univer- 
sality is broken. The spectral statistics in such RMT exhibits unusual 
features like super-strong level repulsion near the origin. 

• Crystallization of eigenvalues 

This phenomenon takes place in the non-perturbative regime k > 1 of 
RMT with log-square confinement and has its counterpart in the "crys- 
tallization" of roots of orthogonal polynomials [18J. So far it does not 
have physical realization, with Calogero- Sutherland model being a clear 
candidate [30] . 

Acknowledgements: I am grateful to my coworkers C.M.Canali, J.T.Chalker, 
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